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I. INTRODUCTION 



£f*l . The problem of detection and measurement optimization at the output of a quantum communication channel, under 
C) its modern interpretation [1], can be described in the following manner. 

Let a quantum channel {H, /5(A)} is given. Here H is the Hubert space where the operators are defined, which 
represent (in accordance with the basic principles of quantum theory) the quantum observables of the channel output; 
/5(A) is the density matrix, which represents the quantum state of the output observables and depends on input 
classical observables A 6 A. The same way as in the classical decision theory [2] a set A of the possible decisions A, 
which are represented now by a final result of some quantum measurement procedure 7. The problem is to find an 
optimal procedure 7 in the standard sense of the decision theory. The latter means that the optimization quality is 
\ represented by a conditional risk function 

r( 7 |A) = f c(\,~X)V 7 (d\\X), (1) 

G\ . J 

' - - 

— , ■ where c(A, A) is the "cost" (pay) function and ■Ry(dAIA) is the conditional probability distribution depending on the 
t-H \ procedure 7 and the conditional channel output state /3(A). 

We suggest the following representation (modified with respect to [1]) of the measurement procedure. In accordance 
with the axioms of quantum theory, a direct quantum measurement of any physical observables represented with a 
commuting set of operators 0, which are defined in the tensor product H ®H' of the given space H and any additional 
auxiliary space H' . The latter introduces some channel- independent quantum system with an arbitrary density matrix 
Qh, p', so that the joint density matrix of the output and the auxiliary subsystems takes the form p(X) <£> p ' . We call this 
indirect measurement. Directly measured variables yield spectral values A £ A, which represents the decisions of 
interest. In accordance with the quantum measurement postulates, we get 



V^{dX\X) =Tr\p(\)®p']E(d\), (2) 
where E is an orthogonal decomposition of unit [3] for operators 

"»h ■ © = / A E(d\). 

A 

Performing in Eq. (2) an averaging over the auxiliary system, we get 

V 1 (d~X\X) = Trp(A)£(dA), (3) 

where 

£(dX) = Tr H 'E{dX)p' (4) 

is the quantum decision function which determines statistical properties of the measurement results. It is the quantum 
analogue of the randomized decision function, that is conditional probability measure p(dX\y), where y are the classical 
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output observables. Here Eq. (3) is analogous to the classical formula V(dX\X) = J p,(dX\y)V (dy\X) . In contrast to p,, 
Eq. (4) defines an operator self-adjoint measure (non-orthogonal, or generalized, decomposition of unit [3]) obeying 
the positivity and normalization conditions: 

£(■) Z 0, J E(dX) = I. (5) 

Thus, optimization over the possible setting (H',p', 0) is reduced to optimization over £ . A similar approach to the 
problems of quantum mathematical statistics is developed in [4], which appeared after the submission of this paper; 
the main results of this paper had been presented in [5] . 

Here only the Bayes type of optimization problems is under study, where an a'priory probability distribution V(dX) 
is known and average risk takes the form 

K = j r(-f\X)V(d\). 

Substituting here expressions (1), (3) and introducing operators of the quantum a'posteriori risk 

R(X) = j c(X,~X)p(X)P(X) (6) 

(which is an analogue of the product R(X\y)V(dy) of the classical a'posteriori risk and the output observables y 
probability distribution), we arrive at the minimization problem 

K(£) = Tr J R(X) £{d~X) = min, (7) 

the solution of which is the optimal decision function £q. 

In this paper the analysis of the possible methods to get the solution of the problem (7) is presented. For gaussian 
channels, for the cases of so called "simple" and quadratic cost functions the optimal decision functions £$ and the 
corresponding optimal measurement procedures are obtained. The most attention is paid to the gaussian density 
matrices 

/5(A) = cxp[r — (x — X) T Q(x - A)], (8) 

where x T = {x\, . . . , x s ) is a given vector of operators in the H space with a non-degenerate (only for simplicity of 
discussion) c-number commutation matrix: 



\X, X 



(xiXj - XjXi) = CI, = 1, ...,s); (9) 



Q > — a given definitely positive symmetric matrix; 

r = lndct 1/2 |2sinh( v /QCv / Q)| (10) 
is the normalization constant. On the base of Eqs. (8), (9) the correlation matrix can be calculated as 

K = (^Tr - Xi)}{ Xj - Xj) + (% - Xj){ Xi - A,)]/5(A)) = ^CcothQC. (11) 

Gaussian density matrices arise at the investigation of the thermodynamically equilibrium linear physical systems 
with a Hamiltonian Ti{X) (that is a total energy represented in terms of coordinates and momenta), which is quadratic 
with coordinates and momenta, i.e. corresponds to a collection of linear oscillators. Typical examples are given by 
vibrational circuits, electromagnetic fields in lines, waveguides and free space. Parameters A in (8) for the most typical 
examples correspond to a non-quantum displacement, due to the presence of a classical signal at their input. Gibbs 
distribution [6] p = exp(r — H/kT) for such systems yields density matrix of the form (8). 
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II. METHOD OF THE SHIFTED A'POSTERIORI RISK 



If the family R of the a'posteriori risk operators is representable in the form 

R(X)=k(X) + A, (12) 
where operator A does not depend on A, then the minimization problem is equivalent to minimization 

Tr^(A) £(d\) = min (13) 

£ 

for the same constraints on £. Such replacement of the minimization problems takes place when the cost function is 
represented in the form 

c(A,A) = c(A,A) + c (A). (14) 

Then the corresponding operators -R(A), A can be expressed via c, c with a help of Eq. (6). 

If the "shifted" family R is commutative, then all operators -R(A) can be represented as functions of some commu- 
tative set of operators y. Then the solution of the problem (13) is easy to obtain in the form 

£ (d\)=6(d\,\ (y)), (15) 
where the scalar function is the indicator function of the subset e: 

S(e,x) = {l f £ % (16) 

Xo(y) — scalar function obtained after minimization with A the matrix elements R yy of the operators -R(A) at the 
^-representation. The corresponding measurement procedure is given then by direct quantum measurement of the 
variables Ao(y), as far as Eq. (15) is an orthogonal expansion of unit for the Xo(y) operators. As far as \o(y) 
are commutative, the decision function (15) can be interpreted as a singular probability distribution of the classical 
type /i(d\\\o) — 5(d\, Ao) and the corresponding measurement as a non-randomized one. Keeping in mind that any 
orthogonal expansion may be written in the form of Eq. (15), we see that any direct measurement is a quantum 
analogue of a non-randomized one. Conversely, the indirect measurement can be called randomized. In the quantum 
case, in contrast to classical Bayes problems, randomization is always necessary when the useful information is 
associated with non-commuting observables and their direct measurement is impossible. 

It is easy to prove, that in the problem of choosing of one of the two possible decisions A = Ai or A = A2 it is ever 

possible to find operator A to fit the demand of commutativity of the shifted risk operators i?(A). Hence, in these 
kind of problems the optimum is ever found on the non-randomized — direct measurements [1]. In all other cases, 
for non-commutative /5(A), it is the most probable that the above method is incapable to provide the commutative 

shifted risk operators R(X). Nevertheless, as it will be shown at Sec. 4, this method may be useful in combination 
with the other methods, which may happen effective in the transformed problem (13). 



III. METHOD OF A LOWER BOUNDARY 

The general sketch of this method is as follows. The minimal average risk 1Z(£) is replaced by some lower boundary 
1Zq(£) ^ 1Z(£). For 1Zq(£) a class X of the decision functions £, which minimize 1Zq(£), is determined: 

TZ (£) = min TZ (£). 

If this class £ contains at least one decision function £ which fits the condition 

K{£ )=K {£ ), 

then £q is the solution of the minimization problem (7), as lZ(£o) ^ Tl{£) for all £. 
One of the simplest lower boundaries is 
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Ti {£)= J r(\)Tr£(\), (17) 

where r(A) are the minimal eigen values of the risks R(X). For the corresponding eigen projectors we have 

P(A)P(A) = r(A)P(A). (18) 

At deduction of Eq. (17) the inequality Tr [R(X) — r(\)I]£(X) ^ is used. For example, if all operators R(X) are 
unitary equivalent all the r(A) = ro are equal. Then with use of Eq. (17) for a finite-dimension space H we get an 
optimal operator measure So in the form of 

£ (d\) = P(X)v{dX), (19) 

if only the positive c-number measure v fits the normalization condition given by the second relation in Eqs. (5). In 
the case of infinite-dimension space H the risk would be infinite, so this case is of no interest. To the latter situation 
the below discussed problem is simplified. 

Let there is a Bayes problem with the simple cost function c(A, A) = — 5(X — A) (A = A are s-dimensional Euclid 
spaces), gaussian family p(X) described by Eqs. (8)— (11) and infinitely wide a'priory distribution p(X) = dV/dX. 
Corresponding a'posteriori risk operators (6) are of the form 

R(X) = -p(A>(A). (20) 
Let us show that in this problem the optimal measure So is given in the form (19) with 

V (dX) = , (21) 

and P(A) are given by the eigen projectors corresponding to to the maximal eigen value po of the density matrices 
/5(A): 

p(A)P(A) = moP(A). (22) 

Now let us divide the space A = A 3 X into regions Aj , which on the one hand have the dimension small enough to 
neglect the dependence of p(X) on A and on the other hand unlimitedly expand with expansion of the p(X) function. 
Then on account of Eq. (20) we can represent risk (7) in the form 

n{S) = -£>(Ai) f Tr p(X)S(dX). (23) 

A, 

For fixed i a'posteriori risks Pj(A) = —p(Xi)p(X) in the right side of Eq. (23) are now unitary equivalent. For the 
subsequent analysis the following asymptotic property is essential for the family of operators /5(A) with A £ A, = 
Aj \ £Aj, that is A belongs to the internal region Aj got by subtraction of a small boundary region <5Aj, the latter 
unlimitedly expanding with p(X) but having a zero relative volume with respect to the one of Aj. For any in advance 
given precision these families may be supposed to act non-trivially only in the subspaces &i(s) C H , e — > 0, beyond 
which, due to relation ||/5(A) , 0|| = ellV'll: it is acceptable to set p{X)ip = 0. That is for unlimited expansion of the 
regions Aj, <5Aj, Aj the subspaces $j are asymptotically orthogonal to each other with any unlimitedly high precision 
given by e. As it is shown in Appendix A, for the discussed family p{X) the orthoprojectors P, onto the introduced 
above subspaces $j can be represented in the form 

P ' = J Pm ^Fi- (24) 

and their dimensionality is iVj = TrPj = det _1 / 2 |27rC|y(Aj) (V — the corresponding phase subspace volume). 

Taking into account equation p(X) = Pj/5(A)p, arising from the above discussion for A e Aj, and neglecting small 
boundary regions SA i: the terms of the sum (23) can be rewritten as 
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TliiSi) = -p(Xi) jTrp(\)£i(d\), (25) 

At 

where £i are the measure £ projection onto $>i fitting the normalization condition 

J £i{d\) = Pi. (26) 

Due to the asymptotic orthogonality of measures £i are independent decompositions of operators P, in $j, so 
minimization of Eq. (23) over £ is equivalent to minimization of Eq. (25) over £%. 

To perform this minimization let us make use of solution (19) with measure v of the form (21) taking into account 
that in accordance with (24) it fits the normalization condition (26). On additional account of independence of measure 
(21) on i this solution (19) defines the optimal measure £o in the total space of A values, so that asymptotically 

Taking into account Eqs. (20) and (22) we get from Eq. (18) 

r(A) = -p(A)/i , (27) 

and for the average risk from (17) or (20), 

K(£ ) = - [ p(X)po 775- = r# • (28) 

J det 1/2 |27rC| det 1/2 |27rC| 

Applying then the zero-temperature method presented in Appendix A it is easy to calculate 

/i = exp|iTr In |2sinhQC| + |QC|] }. (29) 
Correspondingly from Eq. (28) and (29) for C — ► we get the classical limit 1Z = — det 1 ^ 2 (Q/7r). 



IV. METHOD OF REDUCING OF THE A'POSTERIORI RISK OPERATORS TO A FAMILY OF 

UNITARY EQUIVALENT OPERATORS 

In some problems the a'posteriori risks operators may be not unitary equivalent, as it was in Sec. Ill, but with use 
of replacement (12) it is possible to get such a representation for the displaced average risk 1Z(£) = 1Z(£) — Tr A, that 
the method of a lower boundary would be applicable. 

Let us discuss the case of quadratic cost function 

c(A, A) = (A — A) T (A — A) (30) 
(it is no demand to introduce a more general quadratic form, as it is always can be reduced to this one). For this case 

R(X) = J(\-\f(\- X)p(X)V(dX). (31) 

Let us denote 

P = J P {X)V{dX) (32) 
and find such operators u, that Eq. (31) is reduced to 

R(X) =p 1/2 {u-X) T (u- X)p 1/2 + A, (33) 



where A does not depend on A. It is easy to confirm that these operators are the operators of the a'posteriori 
mathematical expectation of variables A: 
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I AF(dA), 



(34) 
where 

F(dX) = f 1/2 p(\) f /2 V{d\) 

is a non-orthogonal decomposition of unit having a meaning of a specific quantum a'posteriori probability distribution 
of variables A . Note, that here introduced analogues of classical a'posteriori characteristics are not unique (e.g. see 
[7]) and the adequate form depends on a specific problem. 

Operator A corresponding to operators (34) in Eq. (33) is equal to 



X T Xp{X)V{dX) -fi /2 u T u~p' 2 i-r,) 
and the corresponding risk is 

1l A =Tr A = Tv(K x -K u ), 

where K\ is a'priori correlation matrix of variables A, for which the equality MA = is suggested; K u denotes the 
condition- less correlation matrix of u depending on the condition- less density matrix (32). 
This way obtained operators of the shifted risk 

R(X) = "p ,2 {u-X) T {u-X)p jL ' 2 

are not unitary equivalent. Although, if in equation TZ(£) = Tr R(X)£(dX) for the average risk we substitute 

~p' 2 E{d~Xyp' 2 = £{d~X), (36) 

then it reduces to 

U(£) = Tr J{u- ~Xf(u - X)£~{d~X), (37) 

where £ is a new positive measure normalized the following way 

£ = p (38) 



instead of the second equality in Eq. (5). For the gaussian density matrix (8) and gaussian V(dX) we have 

p — exp(r — x T Qx) 

with Q = (1/2)C- 1 coth _1 [2(i<: + KxjC- 1 ]. At this from Eq. (34) we get 

u = Vx = K X (K + Kx)- 1 cosh(CQ)x. 

Making use of operators u = K\{K + Kx)^ 1 ^, which determine the optimal operator estimates with no account of 
the commutativity condition [7], we represent this expression in the form 

u = cosh Co Qo uq, 

where C = K X (K + Kx^CiK + Kx)' 1 K x , Qo = [(K + Kx)~ 1 K\\~ 1 Q[K\(K + Kx) -1 ] -1 are, correspondingly the 
commutator of the variables uq and the responding to it Q-matrix in the expression p = cxp(r — u^QoUq) for the 
transformed gaussian density matrix of the output state. 

As a result of the presented transformations, for the multiplier of £ in Eq. (37) we have now unitary equivalent 
operators (u — X) T (ii — A), due to which with use of the lower boundary (17) after replacing £ to £ we immediately 
get the solution in the form of Eq. (19) where projectors P(X) correspond to minimal eigen value ro of operators 
(u-X) T (u-X): 
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(u-X) T (u- A)P(A) = r P(A). (39) 
Corresponding measure f is to be found from the normalization condition (38) which results in equation 



If the equation 



P{~X)u(d~X) = p. (40) 



K u -\\C u \>0 (41) 



holds, then on account of that P(X) is a family of the gaussian density matrices with the correlation matrix Kq = 
(1/2)\C U \ (see Sec. Ill), from Eq. (40) we get measure v in the form of the gaussian probability measure 

v{d\)=det- 1 ' 2 [2^(k u -^\C\ u )] exp[-^F(tf u -^|C u |) _1 A]dA. (42) 

Expressing the original measure £o via the above discussed extreme measure £ o in accordance with Eq. (36) we get 
the optimal quantum decision function 

£ (d~X) = ~f 1/2 P(\yf 1/2 v(d\). (43) 

It is capable to interpret this solution in a simpler equivalent form. For this purpose let us take into account that 
projectors P(X) are of the form tp(X)tp + (X), where are normalized "Glauber" eigen vectors of the operators 

(u- X) T (u-X). Then 

P(X) = [Tr p _1 P(A)] -'f^PixYf 112 (44) 

are also projectors of the form 

P(A) = V {\)<p + (k) = (^+rV)- 1/2 F 1/2 V'). 

- —1/2 

They are, as it is easy to show after multiplication of Eq. (39) both on the left and right by p and applying twice 

„ 1/2 ~ ~ 1/2 ~ ~ 

the commutation rule [7] p u = cxp(—C u Q u ) up (Q u = V~ 1T QV^ 1 ), the eigen projectors of a new generally 
non-Hcrmitian quadratic form: 

[eM-C u Qu)u - A] T [exp(-C„Q u )u - A]P(A) - r P(A). (45) 

Here again, like in Eq. (39), 

r = \Tv\C u \. 

From Eq. (45) it is easy to get that asymptotically with C u — > projectors P(X) and P(X) coincide. With use of 
P(X) Eq. (43) takes the form 

£ (dX) = P(X)D(dX), (46) 

where measure v may be calculated also with use of normalization condition (5) (the second relation). At the limit 
C u — > Eq. (46) differs from optimal measure £o, calculated in Sec. Ill for the case of simple cost function, by the 
change x — > u and applying the corresponding definitions for the projectors P{X) and commutation matrix C. Hence, 
the responding optimal measurement procedures are also identical. 

If condition (41) does not hold then Eq. (42) is not applicable. Not going into detailed analysis of this case, we 
simply mark here, that its qualitative specificity may be revealed by analysis of the transition from the non-degenerate 
matrix K u — (l/2)|C tt | to the degenerate one. At that, over one or several linear combinations of variables u measure 
(42) will be localized at zero (it is the a'priori mathematical expectation). This means that the corresponding variables 
do not need refinement on the base of available information. If matrix K u — (1/2)|C„| is not positive than, from the 
physical qualitative considerations which arc supplied with corresponding mathematical analysis, it is clear that the 
above described degeneration persists. 
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V. METHOD OF INFORMATION CONSTRAINTS 



This method is based [8] on reducing, with use of a small additional term, the linear optimization problem (7) and 
(5) (the second relation) to an appropriate non-linear one which lets using the standard differential variation methods. 
Like in [8] it is possible here to use for this purpose the entropy defined an appropriate way as 



H{£) = J Tr£(dX)\n 



£{dX) 
v(dX) 



(47) 



This expression defines entropy of the relative operator measure £ with respect to the scalar measure v. It generalizes 
usual definition H^l" = j d[i In of the relative entropy of the scalar measure \x with respect to v. Measure v in Eq. 

(47) is to be chosen such a way that the optimal measure £q should be absolutely continuous with respect to v\ to fit 
this condition it is possible to apply the standard measure v{dX) = dX . 

Basing on these considerations, let us switch from the minimization problem (7) to the "regularized" problem 



Tl e {i) = TZ(£)+eH(£) 



mm 

£ 



(48) 



where e > is a small parameter similar to the one used in Appendix A and analogous to a physical temperature; 
1Z(£) is an analogue of the average energy of thermodynamics and H (£ ) is analogue of the entropy. On the total, the 
problem (48) is analogous to the one of the statistical physics involved with the deducing of the Gibbs distribution 
from the minimization of the average energy at the condition of the fixed entropy value. However, this problem differs 
from the thermodynamical one by the operator-valued measure £ and the normalization condition (5) (the second 
relation). The solution £ £ of the problem (48) on account of the normalization condition (5) (the second relation) is 
easy to get with use of the Lagrange multipliers in the form of Gibbs-like distribution 



£ s (dX) = exp 



F e - R(X) 



is(d~X) . 



(49) 



Here an undefined operator F is an operator analogue of a c- number free energy er(e) in a quantum Gibbs distribution 
p = exp[L(e) — H/e}. It is calculated to fit the operator normalization condition 



exp 



F-R(X) 



v(d\) = I. 



(50) 



If the solution F{e) of this equation is acquired the strictly optimal measure £q can be found by calculation of the 
limit £o = lim £ £ . Corresponding to Eq. (49) entropy H = H(£ e ) and average risk 1Z = lZ(£ e ) can be expressed via 

F(e) by formulas which are analogous to the thermodynamical ones: 

dF 

H = Tr — , K = TyF - eH. 

de 

d f 

(To prove this equation one can use the equality — Tr In / £ £ {d\) = 0, which is valid due to the normalization 
condition.) 

At this method the major complication is transferred to the operator equation (50) for the free energy operator F. 
If operator 



p = eln J 



exp 



R{\) 



v(d\) 



does repeatedly commute with all operators R(X), that is 

[[RC\),F },R(X)} = [[R(X),F ],F ]=0 

for all A, then it is a solution of Eq. (50). It is easy to check substituting F into Eq. (50) and using formula 
exp(A + B) = exp(B/2) exp(A) exp(_B/2) for repeatedly commuting operators A and B. Exactly to this situation 
all the problems of Sec. Ill, IV can be reduced. Another possible application of this method may be obtaining 
approximate solutions with use of numeric calculations. 
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VI. PHYSICAL MEASUREMENT PROCEDURE 



Let us confine to the case of quantum decision function 

^>=^>=4^' <5i) 

where P(X) are "vacuum" projectors for operators (x — X) T (x — A) which are quadratic on x: [x,x T ] — CI. As it 
follows from Sec. Ill, IV, this function is optimal for the case of the simple cost function and wide a'priori distribution 
and for the case of quadratic cost function (when x — u) and weak non-commutativity. Let us show that the directly 
measured variables 6 = J XE(dX) have the form 

9 = x + B£, (52) 
where operators £ are copies of x in a space H' a copy of the space H, and B fits the equation 

BCB T = -C. (53) 

The state p' of the additional system is 

p> = P'(0) (54) 

where -P'(O) is a copy of -P(O). 

The proof is given by check of the commutativity of the operators (52) and validness of the Eq. (4). For the 
commutator Cq — [0, T ] on account of Eq. (53) we have C$ = C+BCB T = 0. Then on acount o the commutativity 
of 6 we have 

E(dX) = S(d~X, 6) = 5{dX, x + B£) = (2tt)- s J cxp[^ T (x + s| - X)]dxdX . 

Substituting this expression into Eq. (4), where on account of Eq. (54) is to be represented as Fourier integral 
p' = (2ir)- s f ir(ik)e l ^dk. Then we get 

Tr H >E(dX)p' = d~X(2ir)- 2s JJ n{ix)e l>cT {& -~ X) TY H ,e ll<T B ^+ lkT ^dxdx , 

where the commutativity of x and ^ and possibility of confluence of the exponents of repeatedly commuting operators 
under the Tr operation are taken into account. Applying formula Trexp(?^ T x) = (27r) s det~ 1 ^ 2 \2irC\5(p,) (which is 
easy to prove, e.g. with use of the zero temperature method) and performing integration over k we get, on account 
of the symmetry of ^(ix) and the relation | det B\ = 1 valid due to Eq. (53), we get 

dxdX dX 



Tr H'E(dX)p' = (2 7 r)-Mct- 1 / 2 f n(tB T x)e^ T B ^ d "* X - P{X) 

J dct 27rC 



'det 1/2 |27rC| 

Thus, the presented measurement procedure matches the decision function (51), Q. E. D. 



APPENDIX A: PROOF OF THE SUGGESTED PROPERTIES OF THE FAMILY p(X) FOR AN 

UNLIMITEDLY WIDE VOLUME A 3 A 

Now let us prove the above suggested properties of the family {/5(A), A G A^} for an unlimitedly expanding volumes 
Aj 3 A, SAi. The statements are equal to the relations 

lim(P^ - SijPi) = 0, (Al) 
limp(A) P t = 5 l3 p{X) for A e A r (A2) 

The first one means that Eq. (24) defines projectors onto mutually orthogonal subspaces <E>i, and the second one is 
the condition of the non-trivial action of operators ^(A) only in the subspaces <&i corresponding to Aj 3 X. 
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The proof of the Eqs. (Al) and (A2) for i ^ j is capable to obtain, basing on the explicit expressions of the eigen 
functions Vo(A), which may be used to represent the projectors in the form P(X) = ipo(X)ipQ (X) (this functions are 
frequently called "coherent" or "Glauber" states: for a two-dimensional, or — in other terms — single- mode case 
s = 2 the explicit form of tpo{X) is given, for example, in Ref. [6]). 

Vectors ipo{X) are asymptotically orthogonal to the eigen functions ip n (X') of operators p(X') (n = 0, 1, 2, . . .) if the 
distance A' — A with respect to the norm of commutator C increases to infinity. This is exactly the case of unlimited 
expansion of £Aj. In the case i = j in (Al) and (A2) operators Pj should be represented with the formula (24) 
and, on account of the above discussion, the integration extended to the whole space A. Then Pj will be changed to 
non-orthogonal expansion of unit in the whole space H: 

P(X) . (A3) 

det 1/2 |27rC| 

To get the proof, it is most convenient to use the zero temperature method. For this case it comes to representation 
of projectors P(A) in the form of the limit 

P(A) = limp £ (A) (A4) 

of the density matrix 

(x~X) T (x-X) 



Pe = CXp 



(A5) 



which is a canonical Gibbs distribution for a thermodynamic system with Hamiltonian (x — X) T (x — X) and tem- 
perature e = kT. Let us substitute Eq. (A4) into Eq. (A3) and replace there dX by the gaussian measure 
fi(dX) — exp(— X T K~ 1 X/2)dX, which differs from the probability one by absence of the normalizing multiplier 

det -1 ' 2 (2irK )l ) an falls to dX at the limit — > oo. Then the modified integral is proportional to marginal density 
matrix p(e, K^) (with the proportionality coefficient deX 1 ^ 2 {2-kK ^ / Aet 1 ^ 2 \2i:C\ absent) of the part x of the gaussian 
collection (x, X) with the joint quantum-classical probability distribution /5 e (A)det~ 1 ^ 2 (27rA' ;i )^((iA). Density matrix 
p(e,K^) is gaussian and its calculation algorithm reduces first, in accordance to [7], to calculation of the correlation 
matrix K = K e + K lxi where in accordance to Eq. (11) K £ = (l/2)CcothQC w \C\/2, then calculation from Eq. (11) 
the corresponding matrix 

Q = C- 1 coth- 1 (2KC- 1 ) w (1/2)a7\ 
and with use of Eq. (10) — the corresponding normalization constant 

e f = det 1/2 |2sinh(gC)| w det 1/2 |2 7 rC|/det 1/2 (27rA' AJ ). 
Finally, for the integral of interest we get 



/ 



f,/^ dX ,. ,. dct 1/2 (27rA' M ) /"..,, p(dX) 



det i/z |27rC| if„-ooe-»o det i/J |27rC| J det i/ ^(27rA'^) 



= lim lim t-t- —p(e,K L .)= lim exp — —x K„ x]=I. 

K^ooe^o det 1/2 |27rC| V 2 / 

This proves Eq. (A3) and completes the validation of the above presented solution. 

Submitted 18 I 1972 
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